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Abstract 

We have previously discussed the diffusion limited problem of the bounded one- dimensional mul- 
titrap system where no external fiel is present and pay special attention to the transmission of the 
diffusiny particles through the imperfect traps. We discuss here the case in which an external field 
is included to each trap and find not only the transmission but also the energy associated with the 
diffusing particles in the presence and absence of such a field. From the energy we find the specific 
heat Ch and show that for certain values of the parameters associated with the multitrap system 
it behaves in a manner which is suggestive of phase transition. Moreover, this phase transition is 
demonstrated not only through the conventional single peak at which the specific heat function is 
undifferentiable but also through the less frequent phenomenon of double peaks. 
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I. INTRODUCTION 



The diffusion of classical particles in the presence of traps have been studied by various 



authors 
by others 



The idea of imperfect traps was introduced in [6] and elaborated further 
8| which study, especially, the nearest neighbours distance in the presence of a 
single imperfect trap. We note that the imperfect trap may represent a model for many real 
physical situations that are not easy to tackle directly. Among these one may enumerate all 
the reactions, including the chemical ones, in which the results for the reacting particles can 
not be predicted beforehand. That is, these reactions may result in either the elimination 
of some of the reacting particles or changing the values of some physical parameters that 
are associated with them. Among these parameters one may mention, for example, the 
transmission and absorption coefficients of the diffusing particles, their energy and the other 
variables derived from it such as the specific heat, the entropy etc. In jslQ] the transmission 



properties of these particles in the presence of N imperfect traps, without any external field, 



were intensively discussed for both cases of large and small N. It was found |2( that the 
larger is N the higher is the transmission so that for N — > oo the probability that all the 
particles diffuse through all the traps tends to unity. 

In this work we include an external field to each trap and find the effect of it upon the 
diffusing particles. We note that the case of a bias field in the presence of a single imperfect 
trap at the origin was discussed by Condat et al in [8| which studied among other things the 
effect of a uniform field on the distribution of nearest neighbours distances. In this work we 
introduce the bounded multitrap system and discuss not only the trasmission that results 
from including an external field to each trap but also the energy associated with the diffusing 
particles. This energy is discussed for both cases of the existence and absence of a field and 
we show that it is highly sensitive to the presence of it so that increasing its intensity by a 
small amount may result in a disproportionally large change in the energies of the diffusing 
particles. 

We note that although the physical situation is not strictly of the equilibrium kind we 
follow the tendency of many authors Q| to discuss phase transition and specific heat even at 
situations which are far from equilibrium. In this context we note that it is well known that 
diffusion on a lattice can be mapped onto an equilibrium polymer problem Q|. Note that the 
imperfect traps, which introduces the diffusivity in the multitrap system, are characterized 
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here by a rather small constant trapping rate (denoted by k) as may be seen from the 
numerical work in which we assume k — 1 compared to the value of k = oo for the ideal 
traps. Thus, the bounded multitrap array deviates only slightly from equilibrium. Also, the 
fields added to this equations are entirely uniform as realized from the numerical part of this 
work. We may, therefore, discuss the specific heat and other thermodynamical variables. 

We show in the following that the mentioned large changes in the energy entail, for 
certain values of the parameters associated with the multitrap system, corresponding large 
and even discontinuous changes in the specific heat Ch derived from it. That is, the specific 
heat actually goes through phase transition which is demonstrated not only through the 
conventional undifferentiable single peak but also through double peaks [3, 13 [ill 16, 

Q Q • We note that double peak phase transitions were found in the analogous 
quantum one- dimensional bounded multibarrier potential. Lieung and Neda ^ and also 
Kim et al have found double peaks in the form of the responsive curves which are 
apparently associated with dynamically induced phase transitions. Tanaka et al have 
found such double peaks in antiferromagnetic materials corresponding to magnetic phases 
where the external magnetic field has a corresponding role to the parameter c here. Ko and 
Asakawa have found also double peaks with regards to the phases of the quark-gluon 
plasma which may be thought of as a large number of interactions in a bounded region. 

In Section 2 we represent the appropriate terminology and terms associated with the one- 



dimensional bounded imperfect multitrap system as in |3,Ui|. We then show, using the N 
(2X2) transfer matrices method j^j, 2lj |. that for certain values of the parameters associated 
with the multitrap system the transmission probability thr oug h all the traps tends to unity 
in the presence of an external field. This was shown lloj |. for the absence of field, by 
applying the same method. Using the results of Section 2 we calculate in Section 3 the 
energy associated with the diffusing particles for both cases of the presence and absence 
of a bias field. For this we use the single (4iVX4iV) transfer matrix (as done in j^]) for 
calculating the coefficients of the density of the diffusing particles. In Section 4 we calculate 
the relevant specific heat Ch and show that for certain values of the variables associated with 
the bounded system it behaves in a manner which is suggestive of phase transition. This is 
demonstrated in the presence of an undifferentiable discontinuous peaks in the curves of Ch 
as function of the temperature. Moreover, as noted, we find for some values of the relevant 
variables that this phase transition is demonstrated in the form of double peaks. We then 
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calculate the relevant critical exponents associated with these phase transitions. In Section 
5 we conclude with a brief summary. 



II. THE ONE-DIMENSIONAL BOUNDED IMPERFECT MULTITRAP SYSTEM 

We assume that the imperfect traps, through which the particles diffuse, are all character- 
ized by the same width that depends upon their number N and the length L of the system. 
That is, denoting the total width of the N traps by a and the total interval among them 



re width of each is and the interval between any two neighbours is 
10] the ratio c = - and express a and b in terms of c and the total 



by b we can see that t 
jj. We define as in |9_ 

length L = a + b&sa = , b = . The one- dimensional imperfect multibarrier 
system is assumed to be arrayed along the positive x axis so that it begins from the point 
x = jj = jv^+c) and ends at x = L. We also assume that at t = the particles are con- 
centrated, except for an initial configuration /, at the left hand faces of the N traps. The 
locations of these faces, according to the mentioned arrangement of the multitrap system, 
is given by 

It is assumed that each trap have the same field associated with it so that any particle that 
approach it may be, depending on the direction of the field, either attracted or repelled by it. 
The corresponding initial and boundary balue problem associated with the one-dimensional 
bounded imperfect N trap system is [a 0, Q 



(HyE- X m , t) D pxx{% 1 ^mi t) V p x {x^ X m , t), ^ X, X m ^ L, 

1 < m < N, t>0 
p(x, x m ,0) = e'^d^x -x m ) +f(x,x m ,0), 0<x,x m <L, (2) 



1 < m < N, t>0 
dp V 

Tr\x=x m = (k + —)p(x m ,t), < x m < L, 1 < m < N, t > 0, 

ox D 

where x m denotes the left hand face of the m-th trap (see Eq (JTJ). p(x, x m , t) is the density 
of the diffusing particles and D is the diffusion constant which is considered to have two 
values; Di inside the traps and D a outside them . We note that one may generally 

find in the literature values of D in the range 0.9 > D > 0.3 (see, for example, p. 337 in 
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24j where it is noted that 0.5 is the order of magnitude of D at room temperature and 
atmospheric pressure). We have assigned here for D Q and Di the respective values of 0.8 and 
0.4. p t (x,x m ,t) and p xx (x,x m ,t) are respectively the first time derivative and the second 
spatial derivative. The parameter V is the diffusion velocity that results from the presence 
of the field and its magnitude and sign represent respectively the intensity and direction of 
the field. The second equation of the set (J2J) is the initial condition at t — the first term 
of it signifies, through the S, that there is an initial concentration of the particles at the left 
hand face of the trap located at x — x m , 1 < m < N. The second term at the right hand 
side denotes an additional initial configuration of the density that depends upon x and x m . 
The third equation is the boundary value condition which introduces the velocity V that 
results from the field. The parameter k denotes the degree of imperfection of the traps (in 



it is termed the trapping rate) so that in the limit k — > oo the trap becomes idea 
which case any particle that approaches it is absorbed. We follow the procedure in 
for V = and decompose the set (J2J) into two secondary sets as follows 



in 



PtyXj x mj t) D p xx {x, x m , t) V p x (x, x m , t) , ^ x,x m <J L 

1 < m < N, t>0 



Vx„ 



p(x, x m , 0) = e 2D 5(x — x m ), < x, x m < L, 1 < m < N (3) 

d A - =(k + ^ 

dx Xm v D 



dp V 

(k + — )p(x m ,t), 0<x m <L, l<m<N, t>0 



p t (x, x m , t) = Dp xx (x, x m ,t), < x, x m < L, 1 < m < N, t > 
p{x,x m , 0) = f{x,x, 0), 0<x,x m <L, 1 < m < N (4) 

p(x m , t) = 0, < x rn < L, 1 < m < N, t>0 

The set (jH) is the ideal trap problem in which no external field is present as may be realized 
from the absence of k and V. The set is the imperfect trap problem which includes the 
external field and its initial condition contains only the first term of the corresponding 
condition of the general set (J2J). The solution of (j2J) may be written as in on (see Eq (4) 
there), 

p(x,x m ,t) = Api(x,x m ,t) + Bp 2 {x,x m , t), 0<x,x m <L, (5) 
1 < m < N, t>0 
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where pi(x, x m , t) is the solution of the initial and boundary value problem of the set (J3J) 
and P9.(x, x m , t) is that of (jH). The ideal trap problem of (jlj) is the same as that discussed 
in 3) LLi| so, following the discussion there, we write for p 2 (x,x m ,t) 

Dtix 2 



TXX , 



p 2 (x, x m , t) = sin(— ) exp 



), < x, x m < L, 1 < m < N, 



t > 



(6) 



Regarding pi(x,x m ,t) we may follow Condat et aim [8|, which discuss the problem Q for 
the single trap at the origin, and write the solution for the multitrap 



v V(x-x m ) VH e " ioi 
P!{x, x m ,t) = exp( — ttt)( 



2D 



AD 



yf^Dt 1 2D 



V ^ ■ <>■ ■ - ^^ X -^d + (k + ^)VDi)), (7) 



■ exp(Dt(k + —f + (k + —)(x- x m ))erfc( 



< x, x m < L, 1 < m < X, t > 0, 



where erfc(x) is the complementary error function defined as 25 1 erfc(x) = £°e " 2 c?m. 
In order to adapt the general solution to the bounded X multitrap system we may use, 



as m |y, LLUJ , either the X (2X2) transfer matrices method [20J, |2l| or the equivalent single 
(4XX4X) transfer matrix formalism. We use in this section the X (2X2) transfer matrices 
method for V ^ and write the general matrix equation jj], ^| 



A ™ + i \ (N-l)(a + b) (N-2)(a + b) . 

Bo- l N ] 1 X 



(8) 



>2AT+1 



n(a + b) (n-l)(o + 6) 2 (a + b) a + b I A 1 

" ' 1 X Jl X >••• \ N ' { N ) \ Bi 

where the A's and -B's are respectively the coefficients of the imperfect and ideal trap 
components from Eq (|5)i. The X T's in Eq (jHJ) are all two-dimensional transfer matrices 
that differ from each other by the value of x only jjj. Thus, a representative one which 
relates the two faces of the same trap that are located, for example, at x = x\ and x = x 2 , 
where x 2 > x\, may be written as 



where 



T(xx,x 2 ) 



T 11 (x 1 ,x 2 ) T 12 (xt,x 2 ) 
%i(xi,x 2 ) T 22 (x u x 2 ) 

_ a{D ,xi,t)a(Di,x 2 ,t) 
a(D i ,x 1 ,t)a(D ,x 2 ,t) 



(9) 



(10) 
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T 12 (x 1 ,x 2 ) = (11) 



j- , x x _ rj(Di,x 2 ,i) £(D ,xi,t) _ a(D ,xi,t)£(Di,xi,t) . + 
21 l! 2 f]{D ,x 2 ,t) r){Di,x 1: t) a(D i ,x 1 ,t)r](D i ,x 1 ,t) 
a(D , x 1 ,t) £(A, ^2, *) a(A, ^2, t)£(D , x 2 , t) 



a(Di,X!,t) rj(D ,x 2 ,t) a(D ,x 2 ,t)rj(D ,x 2 ,t) 



(12) 



rr 1 \ v( D o,Xi,t)ri(Di,x 2 ,t) 

T 22 {xi,x 2 ) = -— — — 13 

T](D i ,x 1 ,t)rj(D ,x 2 ,t) 



The parameters a, £, and r\ are given by (compare with the V = case in jjfl, llOj ) 



V(x-x m ) ^..e" 12 ^ V 



a(A x, t) = P2 (x, t) = exp(^ r ^ - i ^)(^ r - (* + ^) • (14) 
• exp( J Dt(A; + ^) 2 + (* + ^)(s - ^))erfc(^^- + (k + ^)VDt)) 



U ' (* + t^) 2 ■ *MDt(k + tttt) 2 + (* + ^)(* - *m)) ■ (15) 



2£>t 7 v 2£r 1 2D' 2D' 



,(Ai., f ] = -- t -iy« (16) 

In order to find the density p(x, x m , t) from Eq (0) at each point x in the multitrap system 
we have to determine the coefficients A and B of the imperfect and ideal trap components 
at this point. If, for example, this point happens to fall at the m-th trap then one have 
to multiply m transfer matrices each of the kind given by Eqs (|1U| ) - (|13|) . We denote the 
components of the two-dimensional matrix that results from such a product by T mn , T mi2 , 
T m21 and T m22 and those of any one of the multiplied m matrices by %^ (I), T 21 (l) and 

T 22 (l) where 1 < I < m. Thus, as realized from Eqs (|10 |) -(|13 |) (and from [10] for V = 0), the 
components T Nll , T Nl2 , T N21 and T N22 may be expressed recursively as 

T Nu = T (JV _ 1)u T u (iV) = ... = T u (N)T n (N - 1) . . . Tn(2)Tn(l) 
Tn 12 = T(jv__i) 12 = . . . = T 2l2 = T 12 (l) = 

T N21 = T (JV „ 1)21 T 22 (iV) + T {N _ lhl T 21 (N) (17) 
T N22 = T (JV „ 1)22 T 22 (iV) = . . . = T 22 (N)T 22 (N - 1) . . . T 22 (2)T 22 (1), 



It is found that the component T/v n tends to unity for growing values of the variables N 
or (and) c. That is, the larger is either the number of traps or the intervals among them 
(or both) the higher is the tendency of the component T/v n to unity. Note that the same 
result have been found also in [9| for V = with respect to the same variables N and c (see 
the discussion there after Eq (30)). We also find regarding the variable V that the lower 
its value is the smaller becomes Tjv n . The same result of small Tjv n is obtained also for 
negative V. We also find that T/v n — > for growing values of the total length L and tends 
to unity for small values of it. Also, Tjv u increases for larger values of either the time t or 

•Em- 

Regarding the component T N22 we find that it does not assume values outside the range 
(0, 1). Also, its factors 7^ 2 satisfy lim7v^oo^22 = for x ~ and lim7v^oo^22 = 1 for x ~ L 
(see in jjj Eq (26) and the unnumbered equation after Eq (27)). Thus, at large N, the 
product T/v 22 also tends to either unity or 0. We also find for small c that the larger is either 
the number of traps N or the time t the more apparently T/v 22 tends to zero and this holds 
also for small values of L. When, however, c increases T/v 22 clearly tends to unity. Note that 
T/v 22 does not depend upon the variable V as seen from Eqs (fT3J) and (fTT)|) . 

The component T^ 21 may generally have any value from (— oo, +oo) but there are specific 
ranges of N, L, c, V, k and t for which T/v 21 tends to zero. Thus, it is found for growing c 
that T/v 21 decreases fastly to zero if the increasing values of c are small and slowly if these 
values are large. For example, increasing c in the range 20 < c < oo causes T/v 21 to decrease 
so slowly that it may be regarded as almost constant over this range. We find that T/v 21 
assumes very high values for either large N or small L but when L grows to a value which 
is comparable to that of iV the component T/v 21 may even decrease to zero. It is also found 
that the larger is V the smaller become the corresponding values of Tjv 21 . For example, for 
V > 24 T/v 21 assumes very small values of the order of 10 -11 but when V decreases below 
20 the value of T N21 becomes very large. Negative values of V have similar effect on T N21 
as the corresponding positive ones. Also, T^ 21 increases for larger values of the time t and 
decreases for small x m . 

We, thus, find that there are ranges of N, L, c, V, k and t for which the components 
of the total matrix Tjv, that result from the iV products at the right hand side of Eq (jHJ), 
assume the values of T/v n = 1, T/v 21 = and T/v 22 = 1. In this case all the diffusing particles 
pass through all the traps in which case the transmission is maximal. That is, the ideal and 
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imperfect trap components of the initial density do not change by the presence of either the 
traps or the fields (or both) in which case the product of the N transfer matrices at the 
right hand side of Eq (JSJ) results in the two-dimensional unity matrix. This is seen in Figure 
1 which is composed of the three Panels A — C which respectively show three-dimensional 
surfaces of the components T/v n , T/v 21 and T/v 22 as functions of c and N. All the three Panels 
are drawn for V = 15, k = 1, D Q = 0.8, Di = 0.4, t = 1 and L = 5 and for the same ranges 
of 15 > iV > 1 and 300 > c > 20. Panel A shows the component T/v n and one may see that 
for large values of c and N T/v n — * 1. The same result is obtained also in Panel C for the 
component T/v 22 whereas in Panel B we see that for large values of c and N the component 
T/v 21 tends to zero. Thus, subtituting these values in Eq (JSJ) and using the fact that one 
always have 0, [l(| T Nl2 = we obtain 

' 1 
1 

That is, all the particles that approach the multitrap system pass through it without any 
decrease in either the ideal or the imperfect trap components of the density. Note that 
for V = we have shown in [l| this unity value of the transmission for the imperfect trap 
component of the density (see Eq (28) there) and the same result is obtained in (lOj] for the 
ideal trap component. 





(18) 



III. THE ENERGY ASSOCIATED WITH THE BOUNDED MULTITRAP SYS- 
TEM 

We, now, discuss the energy associated with the diffusing particles in the presence of 
an external field. We do this by following the conventional discussion one may find in the 
literature regarding diffusive systems in the absence of external fields (see, for example, j3|). 
The presence of the field introduces an additional source of energy (besides that related to 
the diffusion through the traps) which must be taken care of by adding an extra term (to the 
kinetic energy) that depends upon the velocity V. We, thus, assume that the total energy 
is composed of two parts; kinetic and potential where the former results from the diffusive 
motion and the external field and the later from the presence of the traps. In the absence 
of any external field the particles diffuse with an average diffusion velocity vd given, for the 
one-dimensional case, by j3] vd = where D is the diffusion constant. Note that since 
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the two densities inside and outside the traps satisfy Di ^ D Q then also the diffusion velocity 
Vjj and the general density p from Eq (jSJ) satisfy t>£>. ^ vd and p^ 7^ pv . As remarked, the 
natural diffusive motion for V = is towards the positive x axis so when a field is present 
and points in that direction it accelerates the motion of the particles or decelerates it if it 
is oppositely directed. Thus, we may write the kinetic energy of the diffusing particles as 

E Kv¥a (x,x m ,t) = ~p(x,x m ,t)(v 2 D ± V 2 ) = (Ap 1 (x } x m} t) + Bp 2 (x,x m ,t)) ■ 
D V 2 

•(— ± — ), < x, x m < L, 1 < m < N, t > 0, (19) 
t 2i 

where p 2 (a;, x m ,t) and pi(x,x m ,t) are given by Eqs ©-©• V is, as remarked, the velocity 
that results from the external field and the plus and minus signs in front of ^- denote 
respectively that the kinetic energy due to the field is either added for V > or subtracted 
when V < from that due to vd- Note that for large t the kinetic energy from the last 
equation becomes zero since in this case both pi and p 2 vanishes as realized from Eqs ©-(El)- 

By following the conventional discussion of the energy in classical diffusive systems 
we may conclude that the force that acts on the particles is related to the potential energy 
which originates from the presence of the traps. This force is assumed to be proportional 
to the trapping rate k so for very large k (ideal traps) it assumes maximal values and for 
k = (absence of traps) it vanishes. Also, since we have always assigned throughout this 
work a rather small value of unity for k which means that the imperfect traps have weak 
influence upon the particles we correspondingly assume an inverse proportionality of the 
force to the squared distance of the particles from the traps. This means, as remarked, 
that the particles feel the effect of the traps only at small distances from the traps. Thus, 
one may write the force on any particle that results from the trap as F(x,x m ) = — ^-Ar ^ , 
where g is the proportionality constant and the minus sign indicates an attractive force. 
Thus, assuming that each trap serves as a central force source one may find the potential 
energy from 



E P (x, x m ) = - I F(x, x m )dx = - (- - ; )dx = 

J(x r -X m ) J \x r -Xm) \ X X m) 

= -gk{ \ - \ ) = -gk{ 1 ), (20) 

[X X m ) \X r X m ) yX X m ) 

< x, x m < L, 1 < m < N, 
where we assume that the reference point x r is at infinity. Thus, we may write the total 
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energy of the diffusing particles in the presence of field as 



Etotai {V9 , 0) (x, x m , t) = E K(V _, 0) (x, x m , t) + E P (x, x m ) (21) 
< x, x m < L, 1 < m < N, t > 0, 

where E K .(x,x m ,t) and E P (x,x m ) are given respectively by Eqs (fTU|) and (jSDJ). In the 
absence of an external field the total energy is 

1 



E to tai {v=0) (x, x m , t) = E K{V=0) (x, x m , t) + E P (x, x m ) = -p v=0 (x, x m , t)v D + 



+E P (x, x m ) = ^(Ap lv=0 (x,t) + Bp 2 (x,x m ,t)) ■ y, (22) 
< x, x m < L, 1 < m < N, t>0, 

where pi v=0 is the imperfect trap component of the density for V = and is given as 
Pi v=0 = er f(^m) + exp{k 2 Dt + kx) ■ erfc(kVDt + (see Eq (6) in 0). p 2 (x,x m ,t)) 

and Ep(x,x m ) are given respectively by Eqs © and (J2*U|) . The coefficients A and B in 
the last two equations are numerically determined in this section from the single (4NX4N) 
matrix method (see the discussion after Eq (29) in |9j). 

We now show that increasing |V| by even a small amount may change the energy in such 
an unexpected manner that it results, as seen in the following section, in a phase transition 
of the corresponding specific heat C^. This may realized form Panels A-D of Figure 2 which 
all show three dimensional surfaces of the energy E as function of x and c. Panels A and 
B are both drawn for N = 2, x m = j^^y k = t = g = 1, L = 30, D = 0.8, A = 0.4, 
20 > c > 0.5, 40 > x > and differ by the value of V which is 2 for Panel A and 5 for 
B. Thus, by comparing between them one may realize that increasing V by only 3 units 
causes to a disproportionally large increase of E from \E\ w 12 in Panel A to E ~ 14000 
in B. This large jump of the energy entails a corresponding discontinuous change in the 
values of the specific heat Ch which implies, as will be shown, that it goes through a phase 
transition. The same result is obtained also for negative V but compared to V > the 
changes obtained are larger and found at smaller values of negative V. This may be seen 
in Panels C and D of Figure 2 which are both drawn for iV = 5, x m = N 2 ^ c ^ + MxTc) > 
k = g = t = 1, L = 30, D = 0.8, Di = 0.4, 20 > c > 0.5, 40 > x > and differ by 
V which is —0.5 for Panel C and —0.8 for D. Thus, by comparing these two Panels one 
may realize that decreasing V by only three tenths from —0.5 to —0.8 results in a giant 
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SINGLE AND DOUBLE-PEAK PHASE TRANSITION. 
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change of the energy from \E\ m 10 in Panel C to E ~ 10 8 in D. This entails, as will be 
shown in the following section, a corresponding phase transition of the specific heat CV One 
may explain these large changes of the energy by reasoning that increasing the intensity of 
the field \V\ beyond some limit causes the particles to overcome any resistance related to 
the diffusion in the presence of traps. Thus, their energy increases disproportionallv to the 
change of |V| that causes it. A similar behaviour is encountered in Laser tubes (27 1 when 
the pumping energy (field) attains a limit value which causes the intensity of the produced 
light to increase in a phase transitsional manner. 

It is expected, regarding the dependence of the energy upon the trapping rate k, that 
the larger k becomes the more controlled will be the diffusing particles by the traps in 
which case the kinetic (and the total) energy of these particles decrease. This is shown 
from recalling that for N = 2, V = 5, k = 1, t = g = 1, L = 30, D Q = 0.8, 
Di = 0.4, x = n h°^ c \ , 20 > c > 0.5 and 40 > x > we have obtained that the larger 
values of the energy are E ~ 14000. Now, it have been, numerically, found (not shown) 
that if k is raised from k — 1 to k — 5, keeping the values of all the other parameters as 
before, the larger values of the energy are decreased to E « 8000. When k = 10 we find that 
E w 5000 and for k — 15 the larger E further decrease to E m 3000. A similar behaviour of 
decreasing energies for larger k is found also for negative V. 

IV. SINGLE AND DOUBLE-PEAK PHASE TRANSITION IN THE SPECIFIC 
HEAT OF THE BOUNDED MULTITRAP SYSTEM 

The average energy from which one may derive most of the statistical mechanics vari- 
ables such as the specific heat Ch, the free energy F, the entropy S etc may be written as 
< E tota i >= ^ Y°e-^Et ot t a ° tal ) where (3 = kb is the Boltzman constant and T is the tem- 
perature in Kelvin units. Substituting in the former equation the appropriate expression for 
Etotai from Eq (|21|) or (|23J) yields respectively the average energy for the presence or absence 
of an external field. From the expression of the average energy < E tota i > we obtain the 
specific heat Ch 

_ d < E total > _ d EEtotaie-^^ _ 1 2 , 
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Figure 3 shows a three dimensional surface of Ch from Eq (J23J) as function of the ratio c and 
the temperature T where no external field is present in which case the appropriate energy to 
be substituted in the last equation is that from Eq (|22j) . Figure 3 is drawn for N — 2, k — 1, 
L = 30, t = 1, g = 1, x m = jJ^, 20 > c > 0.5 and 10 > T > 0.1. As seen, the height 
of the surface for Ch increases with growing values of the ratio c untill some maximum (not 
shown in the figure). One may also realize that at small values of c and T the surface of 
Ch jumps upward to its local maximal values from which it descends in a similar manner 
to zero. These local maxima are seen to be arrayed along horizontal lines which form, in 
relation to their neighbouring lines, a sharp edge which becomes widened and flattened as c 
grows. 

The epecific heat Ch is certainly undifferentiate at the sharp edge so it goes through 



phase transition at these points 



X 



261 ] . We calculate at the following the critical exponents 



26| associated with this and other discontinuities of CV We find that increasing N or k or 



g, while keeping the values of the other parameters constant, does not cause to any change 
in the form of Ch shown in Figure 3 except to its translation from its position along the c 
axis to a one that tends to be aligned along the T axis. That is, the same surface of Ch is 
rotated in the c — T plane for growing values of N or k or g. If, on the other hand, the value 
of the total length L is simultaneously increased with that of N then the remarked rotation 
of the surface of Ch, obtained for large values of N, is avoided and this surface remain in 
its form and place. Unlike the case of V ^ 0, to be discussed in the following, we find for 
V = that except for the remarked points of discontinuity associated with smaller values 
of c and T there are no other points at which Ch becomes discontinuous. 

We discuss now the specific heat Ch obtained when an external field is included with each 
trap. It is found, as for the V = case, that the corresponding curves of the specific heat 
Ch jumps abruptly from zero for small c and T to their maximal values from which they 
similarly descend to zero. Also, as for Figure 3, these maximal values are arranged along 
lines which form a sharp edge for small c and T which become widened and flattened as c 
increases. The specific heat function is clearly undifferentiable alon g th e sharp edge which 
implies, as for the V = case, that it goes through phase transition [13, 26] at these points. 
But, in contrast to the former case, there exist other points, not at small values of c, for which 
the specific heat goes through phase transition. This is demonstrated in Figure 4 which is 
drawn for V — 5, N — 2, k — g — t — 1, L — 30, x m — jmz^ and which is related not only 
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to the more obvious discontinuity of the two spiky columns at small c and T but also to that 
of the apparently continuous surface at large c's. This is clearly shown in Figure 5 in which 
we isolate from the surface of Figure 4 four curves of the specific heat Ch as function of the 
temerature T for c = 3, 3. 1,3. 2, 3. 3. One can see that each curve of the four shown assumes 
the form of two inverted and indented tooth which are clearly undifferentiable and so they 
constitute a double-peak phase transition. We thus see, as remarked, that the unexpected 
large change in the values of the energy in Panel B of Figure 2, which is drawn for the same 
values of V, N, k, L, g, t, and x m as those of Figures 4-5, is affected through the double 
peak phase transition of the corresponding specific heat of Figure 5. This depends, as noted, 
upon the value of x m = N n+ C ) ' so we ex P ec ^ that changing its value may result in finding 
double peaks phase transitions at other values of c and N. This is indeed the case as we 
find (not shown here), for example, for N = 4, 10.6 > c > 10 and x m = wf^g) + n(i+c) > 
which is the location of the left hand face of the second trap. 

From the former discussion we see that for positive values of V there is associated a single 
peak for the smaller values of c and a double peak for some higher values of it. When we 
consider, however, negative values of V we find that the double peaks generally emerge for the 
smaller values of c. This is demonstrated in the right hand Panel of Figure 6 which shows 2 
curves of the specific heat Ch as function of the temperature T for V = —1.92, g = k = t = 1, 
L = 30, D a = 0.8, Di = 0.4, x m = N ^_ c ^ and for the two values of c = 0.39, 0.4. In this case 
the first peaks of the two curves are small compared to the seconds. The left hand Panel 
shows 8 double peaks curves of the specific heat as function of T for x m = jM^s + jv(i+c) ' 
N = 5, V = -0.37, g = k = t = 1 and c = 1.5 + 0.065 • n, n = 1, 2, ... 7. These phase 
transitions of Ch correspond to the unexpected large change of the energy which is shown in 
Panels C-D of Figure 2 for exactly these values of c, N, t, g, k, x m and in the neighbourhoud 
of V = —0.5. These Panels demonstrate, as remarked, that slightly changing the value of 
V in the neighbourhoud of V — —0.5 by only three tenths changes the larger values of the 
energy from \E\ w 10 1 to E m 10 8 . This change in the energy is demonstrated in the double 
peak phase transition shown in the left hand Panel of Figure 6. Note that all the eight 
first peaks, as well as all the second peaks, touch each other and seem as one curve. The 
appropriate energy to be associated with negative V is that of the expression (|21j) in which 
one should take, as remarked, the minus sign in front of ^r. 

We may suggest an explanation for the occurence of the mentioned large changes in 
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the energy which entail the corresponding discontinuous peaks in the specific heat Ch- We 
confine our attention to the discussed examples of iV = 2 and N = 5 when one respectively 
changes from V = 2 to V = 5 and from V = —0.5 to V — —0.8. As remarked, the change of 
V for N = 2 entails a change in the larger values of the energy from E « 12 to E ~ 14000. 
Looking at the expression (fTT)|) for E one may realize that the large increase in E results 
from a corresponding increase of the imperfect density p\ from Eq (JJJ) (the ideal density p<i 
does not depend upon V and so it does not change with V (see Eq ©))• As seen from Eq 
(J7J) the dependence of p\ upon V is mainly exponential. Thus, when V changes from 2 to 
5 we find for the ratio pi v=5 / Pi v=2 the value of pi v=5 / Pi v=2 = 2.3096 • 10 12 where we use 
the same values used for all the other parameters that lead to Panels A and B of Figure 
2. That is, the density for V = 5 has enormously grows in relation to that for V = 2. 
This is to be compared, for example, to water when one lower its temperature from the 
gaseous state to the liquid one in which case the density of the water molecules grows in a 
phase transitional manner. This behaviour is repeated when one continues to decrease the 
temperature to 0°C from the liquid state to the solid one in which case the density of the 
water molecules increases again in a phase transition manner. We have mentioned in the 
previous section the example of laser tubes for which the intensity of light increases greatly 
when the pumping energy (corresponding to V here) attains a specific value. This occurs 
because a macroscopic aggregate of atoms have been transferred by the pumping energy into 
the appropriate laser state 11]. A similar grow of the density occurs, as remarked, also here 
when one increases for N = 2 the velocity V from 2 to 5. 

The double peaks shown in Figure 6 which are associated with small negative V may 
be explained by noting that the external diffusion constant employed here is D Q is 0.8 (the 
internal diffusion constant is even lower Di = 0.4). That is, when one turns on an external 
field which is directed opposite to the diffusive motion then when the value of this field 
becomes V = —0.8 it actually neutralizes and cancels the influence of the traps on the 
particles so that their energy becomes very large. 

The appearance of the double peaks for these values of V demonstrates further (more 
than the single peak) the large change that the density and the energy have passed through 
when V changes as described. We must, however, note that these peaks depend not only 
upon V but also upon the other parameters, such as N, x m , t etc, that control the behaviour 
of E. 
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one may calculate the related critical exponents x [26| associated with these phase tran- 
sitions by using the following equation in the neighbourhood of the critical temperature T c 

Q 

C h (e) = A + Be x , (24) 

where e = r ~ Tc and A, B are constant. The first order derivative of C/j(e) diverges at T = T c 
so the critical exponent x ma y be obtained from [26j 

ln\C h (e)\ ln\§\ 
X = 1 + hm = 1 + hm / (25) 

where the unity term denotes the first order derivative Cfc(e) of the specific heat with respect 
to e. The value of \ may be obtained by plotting the curve of Cvj(e) in the close neighbour- 
hood of T = T c and one can see from Figures 3-6 that \ assumes different values. Thus, 
assigning to the constants A, B the respective values of and 1 and plotting, as remarked, 
the graphs of C^(e) in the immediate neighbourhoud of the single peaks in Figure 3, which 
are located at small c and T, one may calculate Xi using Eqs pij) -(|2*5 ]) . as x ~ §■ This value 
changes with respect to the single peaks of Figure 4 which are also located at small c and 
T. That is, repeating the same procedure one may obtain the value of x ~ §• F° r the eight 
first and second peaks in the left hand Panel of Figure 6 we find the respective values of 
X ~ 3 and x ~ §' ^ n e x °f the double peaks of Figure 5 and also of the first peaks in the 
right hand Panel of Figure 6 is x ~ §• 

Comparing the phase transition behaviour of the bounded one- dimensional multitrap 
system to that of the corresponding quantum array of the bounded one-dimensional multi- 
barrier potential ^| one may notice the following similarities and differences; The specific 
heats of both systems exhibit the same discontinuous jump at small values of T and c but 
whereas in the multitrap system Ch decreases, for growing values of T, to zero the corre- 
sponding quantum Ch does not vanish but tends, for large T, to a finite value (see Figures 



1-7 in 



14|). Also, the phenomenon of the double peaks at which the specific heat Ch is 



undifferentiable are discernable in both systems. 

The variation of the critical exponent x may be explained by noticing from Figures 3-6 
that the different peaks shown in these figures correspond to different values of N and c. 
Remembering that N and c respectively denote the number of traps and the ratio of their 
total interval to their total width one may realize that they, actually, control the shape 
and form of the multitrap system. That is, Figures 3-6 with the different values of N 
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and c actually correspond to different systems through which the particles pass and not 
to different parameters of the same system. Thus, one may expect different values of the 
critical exponent to be associated with these different systems. We note, however, that the 
difference between these values is not large. 

Using the expressions (j21j) -()22 j) for the energy one may obtain the other variables of 
statistical mechanics. For example, the free energy F is calculated, for either the presence 
or absence of the external field, from 2^] F = —khTlnjJ^, e~^ Etotal ). Using the last equation 
one may write the entropy S for the multitrap system 



24] 



S=~ = ^Tln($>-^-)) = h\n(£e-^) + h ^^%^ (26) 

As for the specific heat Ch one may draw S for different values of the parameters N, k, g, 
x m , t, L and c. If, for example, we draw the surface of S for exactly the same values of the 
mentioned parameters as those of Figure 4 one may see (not shown here) two separate lobes 
for small c and T which correspond to the two spiky columns of Figure 4. 

V. CONCLUDING REMARKS 



We have discussed the diffusion limited problem related to the bounded one-dimensional 
imperfect multitrap system in the presence of external field V. The analytical methods 
previously used 0, [l0| to discuss the transmission of the particles through this system in 
the absence of an external field were also used here in the presence of it. Thus, the N (2X2) 
transfer matrices were used, as in 19|, for discussing the transmission through the multitrap 
system for V ^ and the single (4NX4N) transfer matrix for studying the energy and 
the corresponding specific heat CV It has been found, as for the V = case in |9j, that 
or certain values of the parameters associated with the system the transmission coefficient 
i3 Q| of the diffusing particles tends to unity when N and c become large in which case 
all the particles diffuse through all the traps. This has been shown not only for positive V 
which pushes the particles towards the traps but also for negative values of it that repell 
the particles towards the negative direction of the x axis. 

The unique characteristics of the mutitrap system become more unexpected regarding 
the energy of the diffusing particles in the presence of an external field. Thus, it has been 
found that increasing V, for either positive or negative values of it, by even a small amount 
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results in a disproportionally large increase in the energy of the diffusing particles so that 
trying to calculate the related specific heat Ch we find that it goes through phase transition. 
Moreover, for certain values of the parameters associated with the multitrap system, such as 
its total length L, the number of traps N, the ratio c, the time t, the location x m at which 
the particles are initially concentrated and the field V, one may find that the mentioned 
phase transition is demonstrated in the form of a double peak. The value of the related 
critical exponents associated with these phase transitions were found to vary between | and 

3 
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FIG. 1: The three Panels A — C show three-dimensional surfaces of the components T/v n , Tjv 21 
and T/v 22 from Eqs (|17|1 as function of the ratio c and the number of traps N in the ranges of 
300 > c > 20 and 15 > N > 1. The other variables are assigned the following values; D Q = 
0.8, Di = 0.4, t = k = 1, L = 5. The parameters c and iV are obviously dimensionless and one may 
realize from Eqs (|5l)- lfT7|) that the three components Tjv 115 T/v 21 and T/v 22 are also dimensionless. 
As seen, the components Tjv n and Tjv 22 tend to unity for the larger values of c and N and the 
component Tjy 21 tends to zero. This signifies that all the diffusing particles pass the multitrap 
system (see Eq l(TS|) ). 
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FIG. 2: Each of the four Panels shows a three-dimensional surface of the energy from Eq (|21[l as 
function of cc and the ratio c. The length x is given in units of cm, c is dimensionless and the 
energy E is given in units of ergs. Panels A — B are both drawn for N = 2, D a = 0.8, D{ = 0.4, 
t = k = g = 1, L = 30, and x m = N ^ c +C ^ but V = 2 for Panel A and V = 5 for B. Note that 
by increasing the value of V from 2 to 5 results in an unexpected large change of the energy from 
\E\ ~ 10 to E ~ 14000. Similar results are shown in Panels C — D which are both drawn for 
g = k = t = 1, L = 30, D = 0.8, A = 0.4, N = 5 and x m = jm^\ + N (l+c) but V = ~ - 5 for 
Panel C and V = —0.8 for D. Note the giant change from \E\ ~ 18 in Panel C to E ~ 5 • 10 8 in 
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FIG. 3: A three-dimensional surface of the specific heat Ch from Eq Q23JI as function of the ratio 
c and the temerature T and in the absence of an external field. The units of Ch and T in this 
figure and in Figures 4-6 are [ d ^ ee ] and Kelvin degree respectively and c is dimensionless. The 
appropriate expression for the energy substituted in Eq (|23|) is that from Eq (|22|) . The figure is 
drawn for N = 2, D a = 0.8, Di = 0.4, L = 30, k = g = t = 1 and x m = N ^_ c ^ ■ Note the sharp 
edge of the surface for small c and T which becomes widened and flattened as c increases. 
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FIG. 4: The specific heat Ch from Eq ()23|) as function of c and T. This Figure is drawn for V = 5, 
N = 2, D Q = 0.8, Di = 0.4, L = 30, k = g = t = 1 and x m = N ^ +c ^ ■ Note that these are the values 
for which Panel B of Figure 2 is drawn. From that Panel we see that changing the value of V fron 
2 to 5 results in a large change of the energy. This change is demonstrated in the discontinuity of 
Ch for either small c as in here or for larger values of it as in Figure 5. The units of Ch and T are, 
as remarked, [ d ^ ee ] and Kelvin degree respectively and c is dimensionless. 
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FIG. 5: The figure shows four different curves of the specific heat Ch in units of [ d ^ ee ] as function 
of the temperature T in units of Kelvin for N = 2, V = 5, D Q = 0.8, Dj = 0.4, g = k = t = 1, 
L = 30, x m = jv(i+ e ) and for the following 4 values of c; c = 3, 3.1, 3.2, 3.3. The double peaks are 
clearly seen in each curve. Note that the spiky forms of Figure 4 are obtained for exactly the same 
values as in this figure except that c is smaller. Thus, one may conclude from Figures 4-5 that for 
positive V there exist single peaks for small c and double peaks for the larger values of it. 
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FIG. 6: The two Panels show double peaks of the specific heat Ch as function of T for negative 

V and D = 0.8, A = 0.4, g = k = t = 1, and L = 30. The Panel at the left shows eight 
double-peak curves for N = 5, V = —0.37, x m = j^+cj + N(i+c) an< ^ ^ values of c = 
1.5 + 0.065 ■ n, n = 1,2, ... 7. The Panel at the right shows two double-peak curves for N = 4, 

V = —1.92, x m = N ({ c +C ^ and the two values of c = 0.39, 0.4. C\ and T are given in units of [ d ^ ee ] 
and kelvin degree respectively. 



